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We study dynamic fluctuations in non-disordered finite dimensional ferromagnetic systems 
quenched to the critical point and the low-temperature phase. We investigate the fluctuations 
of two two-time quantities, called x an d C, the averages of which yield the self linear response and 
correlation function. We introduce a restricted average of the x' s i summing over all configurations 
with a given value of C. We find that the restricted average {x)c obeys a scaling form, and that 
the slope of the scaling function approaches the universal value Xoo of the limiting effective tem- 
perature in the long-time limit and for C — > 0. Our results tend to confirm the expectation that 
time-reparametrization invariance is not realized in coarsening systems at criticality. Finally, we 
discuss possible experimental tests of our proposal. 

PACS: 05.70.Ln, 75.40.Gb, 05.40.-a 

I. INTRODUCTION 

Fluctuation-dissipation relations (FDRs), namely model- independent relations between linear response functions 
and correlation functions, have been extensively investigated in systems that relax slowly out of equilibrium [l|, 0, 0, 
EL IS IE 0- Special emphasis was set on the analysis of aging cases. For spin systems the impulsive auto-response 
function, describing the effect of a perturbing magnetic field hi(t') acting on site i at time t 1 on the magnetization 
(si(t)) on the same site i at the later time t > t' , is 

The (spatially averaged) integrated auto-response function, or dynamic susceptibility, is 

(x)(M»)4E t {Ri)(t,t')dt' , (2) 

where N is the number of spins in the system. For ferromagnetic coarsening dynamics it is convenient to consider a 
perturbation that is not correlated with the equilibrium ordered states and quenched random fields hi are typically 
used. For instance, in the case in which a bimodal random field, hi = hei, with h the amplitude and e$ = ±1 with 
probability a half is applied, the susceptibility J2]) can be cast as 

(X) (t, t w ) = hm - ^ , (3) 

i 

where s 1 ^ (t) is the value that the i-th spin takes at time t in a trajectory in which the perturbation was switched 
on from t w onwards, and Si(t) is the value of the spin in a freely evolving trajectory. Other choices of random 
perturbations with, for instance, finite spatial correlation are also of interest [8( but we do not use them here. The 
average (. . .) is taken over all possible initial conditions, thermal histories, the random field and quenched disorder (if 
present). With this notation, (x) and (R), we make explicit the fact that we average over all sources of fluctuations. 
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In equilibrium, (x) and the auto-correlation function 

(C)(t,t w ) = i^( Si (t)si(t„)) - -^£>i(*)>^X>i(*«,)> , (4) 

computed with the same complete average, depend only upon the difference t— t w , due to stationarity, and they are 
related through the fluctuation dissipation theorem (FDT): T(x)(t - t w ) = {C)(t,t) - (C)(t - t w ) = 1 - {C)(t - t w ) 
(we consider here and in the following unitary modulus spins). In generic non-equilibrium states, (C) is no longer 
stationary but it is, usually, a monotonically decaying function of t. One can then invert the relation between t and 
(C) to obtain 

{x){%t w )=x{{C),t w ) ■ (5) 
Also quite generally x((C)?^«>) looses the t w dependence at large t w and the limiting form 

2«C»= lim x({C),t w ) (6) 

t w — >oc 

is a non-trivial function of (C) l7j. Moreover, the slope X((C)) = —Tdx((C))/d(C) allows for the definition of an 
effective temperature |6| through T e //((C)) = T/X{{C)). For coarsening systems quenched to the critical point the 
limiting value 

Xoc = -T lim lim ^^77^- (7) 

is of particular relevance [2, H, 0, H, H EH EL E2] ■ Note that taking (C) — > first implies that one takes t — > 00 before 
t w — ► 00. 

What discussed insofar shows that important properties of the system, the effective temperature in particular, are 
encoded in the dependence of (x) on (C), with t and t w being parameters the variation of which merely allows one to 
scan sectors with different values of the relevant quantity (C). In this paper, we elaborate on this idea by studying 
the fluctuations of these two-time functions. In complete generality, for a given initial condition, thermal noise and 
random field realization we consider the fluctuating quantities 



I N 

C(t,t w ) = -^ Si (t) Si (t w ) , (8) 



i=i 



i=l 

without any averaging (we dropped the second term on the right hand side of Eq. ^ since it is negligible for large N 
in the cases with (s^ = considered in the following). The joint probability distribution P(C, x) has been studied 
in disordered spin models [HI, [13, some kinetically constrained spin systems fl5^ . and the 0{M) ferromagnetic 
coarsening in the infinite N limit jld| and including 1/M corrections [TtJ with the aim of checking predictions from 
the time-reparametrization invariance scenario of glassy dynamics [l8| . 

The analysis of the statistics of two-time fluctuations can be more conveniently carried over if some fluctuations 
are damped by performing a partial averaging of x, in the following referred to as restricted average. Specifically, we 
introduce the quantity 

M c (u.)^falE (( ' ; ' ( "7" )K>c M 

i 

where the global average (. . .) of Eq. ^ is replaced by a restricted average (. . .)c over configurations with a given 
value of C as sketched in Fig. [T] As will be explained in the following, considering restricted averages greatly simplifies 
the analysis still retaining some basic physical information on the relation between x and C . Since, even for t and t w 
fixed, C is a fluctuating quantity, the restricted averaging procedure allows us to explore the behavior of (x)c as a 
function of C, similarly to what is done with the globally averaged quantities (x) and (C) by varying t and t w (notice 
that, differently from globally averaged quantities, the plot of (x)c vs C has also a negative branch). The physical 
idea inspiring this analysis is that the dependence of (x)c on C should bear the same information as x((C)) on some 
properties of the system, in particular on the limiting value T/X 00l the effective temperature. 
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Let us state the problem more precisely. Generally (x)c depends on C, t and t w (see Fig. Q}. Taking advantage 
of the monotonic properties of (C) and (x) as functions of times, one can replace the temporal dependencies in favor 
of (C) and (%). The restricted averaged quantity (x)c can then be recast as (x)c = Xc{C, (C), (x))- In the limit 
t w — > oo, the dependence on (x) can be dropped since in the models considered (x) approaches x((C)), see Eq. 
and it becomes redundant. Hence 

lim (x)c(t, M - , lim (x)c(C, (C)) . (11) 

Now, it is convenient to extract a factor x((C)) ~ the integrated fully averaged linear response- from the right-hand- 
side, and to use the natural variable C — (C) describing the fluctuations of C around the average to write: 



lim ( X )c(t,t w ) = x((C)) (x)c(C-(C),(C)) (12) 



[for simplicity we use the same symbols for the two different functions (x)c m Eqs. (fTTj) and l|12p]. Equation (fT2"|) is 
simply a rewriting of (x)c{titw) m terms of the more natural variables C and (C). From this point on we proceed 
by using some physically motivated assumptions, the soundness of which will be tested in Sec. IIII1 Let us notice first 
that in Eq. (|12[) we are left with a two-parameter dependence in (x)cj while x((C)) depends on a single parameter. 
The guiding idea of an analogy between (x)c(M«0 and x((C)) discussed above suggests that C and (C) may enter 
(x)c iu a particular combination, thus reducing the number of parameters to one, as for x((C))- Since C has the 
upper bound C(t,t) = 1, the natural scale of fluctuations is 1 — (C). Therefore, we make the following scaling Ansatz 

lim ( X ) c (t, t w ) = x«C» / (tf^gf) , (13) 



C-oo W " V ' J \1-{C) 

where f(x) is a scaling function from which, we conjecture, one can extract Xqo in the region C = 0. In this paper 
we check this conjecture in the Ising model in d = 1,2,3 quenched to the critical temperature or below. We find 
that the scaling form (fTB")) is verified, and that the slope of f(x) evaluated at xq = —(C) /(I — (C)), the x value that 
corresponds to C — 0, yields the limiting X^: 

-^|«o=*co. (14) 

This claim is done asymptotically and we discuss its implications and how to take this limit in the body of the paper. 

This Article is organized as follows. In Sec. [II] we overview what is known about the scaling behavior of (x) and 
(C) in quenched ferromagnetic systems. In Sec. IIII1 after defining the restricted average and the methods used to 
measure it, we check the validity of the Ansatz (IT5]) and (fl"i|) in various systems. Specifically, in Sec. lIII Al we consider 
the Ising model in d = 1 quenched to T = in an independent interface approximation (Sec. IIII A 1[) and by means of 
numerical simulations (Sec. IIII A 2|) . In Sec. IIII B I we then study the Ising model in d = 2, 3 quenched to T c with an 
analytical Gaussian approximation (Sec. IIII B lj) and numerically (Sec. IIII B 2[) . Finally, the case of a quench of the 
d > 1 Ising model below T c is considered in Sec. lIII Cl We discuss the results, the relations with the reparametrization 
invariance symmetry and some open problems in Sec. IIVI 



II. THE AVERAGED LINEAR RESPONSE AND CORRELATION IN NON-DISORDERED 

COARSENING SYSTEMS 

Before entering the field of fluctuations, it is useful to overview the pattern of scaling behavior of the fully averaged 
quantities (x) and (C) which are quite well understood in the relatively simple case of clean (without quenched 
disorder) ferromagnetic systems. 

In the case of quenches to the critical point of a scalar ferromagnetic model the averaged two-time functions satisfy 
the scaling forms 

(C)(t,t w ) = t- a g(t/t w ,t Q /t w ) , (15) 
{x)(t,t w ) = x eq -t~ a g(t/t w ,t /t w ) , (16) 

where x eq = /? is the static equilibrium susceptibility, a — 0.115 in d = 2 [li| [2(| and a = 0.506 in d = 3 [2(|. In 
Eqs. (fT5) and (fTo) t is a microscopic time needed to regularize (C) and (x) at t/t w — 1 ensuring (C)(t,t) — 1 and 
(x)(£, t) = 0. In the asymptotic limit t w — > oo for any t/t w fixed the averaged correlation vanishes and the integrated 
linear response reaches the equilibrium value x eq = P due to the t~ a prefactors. Equations (fT5)) and fTo]) imply @ 



TX«C)) = 1 - (C) , 



(17) 
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FIG. 1: The procedure used to compute (x)c(t,tw) is illustrated for the Id Ising model quenched to T = 0. The dashed line 
is the curve x((C)) an d the bold symbols are the points of coordinates ((C), (x)) for the two cases considered. For a given 
choice of t w and t, that in the figure are t m — 10 and t — 20 (right) or t — 100 (left), in units of Monte Carlo steps (MCs), one 
computes the joint probability distribution P(C,x)- Some contour levels are sketched as ellipsoidal curves in the figure. The 
average of x at fixed C (along the vertical direction) is then calculated. By varying C one obtains the solid curves depicted in 
the figure. In the case t = 20 (right) we show the restricted averaged integrated response obtained using Eq. (|10|) or by means 
of the field- free method based on Eq. (|23[1 . as discussed in the text. The two curves superpose and are almost indistinguishable. 



and for the explicit forms of g and g found, interestingly, the limiting slope Xoo takes a non-trivial universal value 
Xqo < 1, that if interpreted as yielding an inverse effective temperature, T/Xoc, signals the presence of a higher 
effective temperature than the bath one in the peculiar limit of Eq. Q. The value of has been determined using 
field theoretical techniques up to second order in e = 4 — d [3]; one finds — 0.429(6) and Xoc = 0.30(5) in d = 3 
and d = 2, respectively. Numerical calculations yield X^ ~ 0.4 in d = 3 pi and X^ ~ 0.33 — 0.34 in d = 2 [2lj . 

For quenches below the critical point, two time quantities split up [22| into a stationary (quasi-equilibrium) part 
and an aging part 

(C)(t,t w ) = (C) st (t-t w ) + (C) as (t,t w ) , (18) 
(X)(t,t w ) = (x) st (t-t w ) + ( X ) a9 (t,t w ), (19) 

with (x) st and (C) st related by the equilibrium FDT. This means that (x) st saturates to the static susceptibility 
X eq — /3(1 — M 2 ), M being the equilibrium magnetization density below T c , in the finite characteristic time of the 
equilibrium state. The aging parts obey the scaling forms [3, [23[ 

(C) ag (t,t w ) = g(t/t w ) , (20) 

(x) ag (t,t w ) = tZ, a W/t w ) , (21) 

for large t w . Above the lower critical dimension d^ one has a > 0, implying lim tm _ +00 (x)(i, t w ) = (x) st (t — tw) and 

1X{{^)) - | !_ M 2 for ( C )<M 2 . {22 > 

Therefore X^ — 0. Let us stress that x((C)) is determined by the stationary part of the response function only, with 
the aging parts producing finite time corrections which depend on (C) and t w . 

The same scaling structure (|T8]) - (f^T|) holds in the case of systems at the lower critical dimension quenched to T = 0, 
with the difference that the stationary parts vanish [2~i| and a — 0. This implies that y((C)) is, in this case, a property 
of the aging regime. Its exact computation yields a non-trivial form with X^ — 1/2 [251 ]. 
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III. RESULTS FOR ( x )c 

In this section we study the behavior of the restricted averaged response (x)c hi different systems. 

Reaching the vanishing applied held limit in the calculation of the response function is a difficult task that has 
been discussed in numerous studies of the fully averaged response. A variety of methods that avoid applyin g a field 
and transform the globally averaged linear response into correlation functions have been proposed and tested |26l . H3] • 
Elaborating on these ideas, we argue that not only (x), but also the restricted average, (x)c> can be computed over 
unperturbed trajectories. Actually, following the same line of reasoning exposed in [27T ] but taking at the end of the 
calculation just the restricted average over the x's leads to 

jv t 

2T( X )c(t, t w )=l- C(t, t w ) + - ]T Bi(t')) c , (23) 

where 

Bi = 2s lWl , (24) 

Wi is the transition rate for flipping the spin Si and t' runs over elementary moves. Notice that this relation between 
the response function and unperturbed quantities is only valid for the averages (global or restricted) of %, while the 
relation between the right hand side of Eq. (|23p and the fully fluctuating quantity \ is n °tj i n principle, known. 
The availability of the fluctuation-dissipation relation (|23"|) , therefore, is one of the great advantages of dealing with 
restricted averages instead of considering directly x- 

In numerical calculations one evolves many different realizations of the system up to time t grouping them according 
to the value of the fluctuating global overlap C . (x)c can then be computed in two equivalent ways: 

i) At time t w a replica of the system is created on which the perturbation is switched on, and the fluctuating 
quantity x is computed through Eq. ((3j. With the set of points (C, x) °ne then computes the restricted average (x)c 
by averaging the x's over all instances with the same C, as described in Fig. [T] A dependence on the magnitude of 
the applied field remains and one is interested in the h — > limit. 

ii) The right hand side of Eq. (|23|) (which is itself a fluctuating quantity) is computed on the unperturbed trajectory 
and hence averaged over realizations with the same C. There is no applied field in this case. 

As shown in Fig. [T]the two methods give identical results, but the second one is much more efficient computationally 
with, moreover, the built-in limit h — > 0. In the following, therefore, we shall compute (x)c using Eq. (|23p . 

In the definitions above the C and x are summed over all spins in the sample and this poses a problem. The use of 
large system sizes £, N ^> 1, is needed to avoid too important finite-size effects. But for very large systems significant 
fluctuations of C are rare and one can only access tiny variations around (C). In order to avoid this difficulty we prefer 
to collect the statistics over subsystems with, say, linear size £, each containing n — £ d spins, and then compute the 
restricted average through Eq. (|23|) with n in place of N. Ideally, the coarse graining length £ should be much larger 
than the lattice spacing t 3> S. Actually, the results of the analysis carried out in this paper become independent of 
the coarse-graining length for large £, but if i is too small hnite-size effects affect the statistics of C and (x)c much 
in the same way as a small C would introduce corrections with respect to the thermodynamic limit C — > oo in the 
global quantities. Then, a convenient choice of the coarse graining length is to fix I to be sufficiently large so as to 
avoid significant finite size effects, but not too large either, otherwise fluctuations would be too rare. Since finite 
size effects in coarsening systems occur on length-scales of the order of the growing length L(t), a convenient and 
realizable choice of t is 

S < L(t) < £ < C. (25) 

A similar choice was made in [3] and [28| for the study of fluctuations in the 3d Edwards- Anderson spin-glass and 
the random field Ising model, respectively The scaling of the correlation fluctuations with the additional variable 
£/L{t) during coarsening was also discussed in (28|. In this paper we are not interested in checking this kind of scaling 
but we simply use £ as a magnifying glass to tune the extent of fluctuations under investigation. 



A. d = 1 Ising model quenched to T = 

We here focus on the one dimensional Ising model prepared in an initially disordered configuration at infinite 
temperature and evolving with Glauber dynamics at T = 0. The dynamics is then simply Brownian diffusion of 
interfaces with annihilation upon meeting. The averaged aging properties of this model are well established [9(. In 
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two recent papers [29| the study of multi-point correlation functions revealed the existence of dynamic heterogeneities. 
Here, we extend these studies to analyze the relation between the fluctuations of two-time quantities that once averaged 
become the linear response and correlation. 



1. Non interacting interfaces approximation 

Here we derive the expression for (x)c i n an approximation in which interface annihilation is neglected but interfaces 
travel a long distance in the interval t — t w . This is expected to be correct in the regime t w 3> to (interfaces are 
very distant in the sample) and t/t w > 1 (the chosen interface travels a long distance in the t — t w interval and C is 
significantly different from 1 although the limit t/t w — > oo is not reached). Interestingly enough, we shall see that the 
results derived in this limit capture the basic features of the fully interacting system for all t/t w , even for t/t w 1 
when interactions between interfaces should become important. 

For a single interface one can compute (x)c exactly (see Appendix I for details). The result is 

2T(x)c = l-C+(D) c , (26) 



where 

(D)c = ^J2 X! (2n e + l+v w ) P c (n c ,v w ) . (27) 



2 



N 

n c — v w — ±1 



n c is the number of times the interface passes, in the interval through its final position at time t, i.e. x(t) 

(see Fig. [2]). v w — +1 (v w = —1) if the first move of the interface (at time t w ) is towards (away from) the final 
position, x(t). In the sketch in Fig. [5] this means that v w = 1 (v w = —1) if the upper configuration moves to the 
right (left) in the next time-step. Pc(n Cl v w ) is the probability of finding a particular realization of n c and v w in the 
restricted ensemble with a fixed value of C. Notice that 1 — C = 2Ax/N is entirely determined by the initial and final 
position of the interface in terms of the distance Ax = x(t) — x(t w ) traveled. A dependence on x(t) and x(t w ) also 
occurs in (D)c through v w , because the probability of moving toward the final position is always larger than that of 
moving in the opposite direction. However, for sufficiently large values of 1 — C, when the interface has traveled a long 
distance, the sum is completely dominated by the term Y^ n =o (2n c + l)Pc( n c), where Pc(n c ) = J2 V =±i Pcfaa v w)- 
Moreover, in the same limit, Pc(n c ) itself does not depend on the restriction C, because the position of the interface 
becomes independent on the number of times the interfaces passes through x(t). In this limit, therefore, (D)c does 
not dependent on x(t) and x{t w ) and, hence, it is independent of C. (D)c can be determined by using the 'sum rule' 
J2c(x)cP(C) = (x), where P(C) is the probability distribution of C. This yields (D) c = 2T( X ) - 1 + (C). Hence 
one has T(x)c — T(x) — (1/2)(C — (C)). Recalling [30( that for a single interface T{x) = 1 — C, (x)c can be cast in 
the scaling form (fT5)) 

nmJ X )c(t,t w )=Z(C))f(^^\ with f{x) = l- l -x. (28) 

The linear dependence of the function / on its argument implies that its slope is always equal to —1/2 = X^, which 
agrees with Eq. (|14[) . This result, however, is expected to apply only within the limits of the approximation, that is 
to say, when C is not too close to 1 nor to —1. We put this prediction to the numerical test in the next subsection. 



2. Numerical results for the full model 



In this section we study numerically the behavior of the fully interacting model. Here, and in the following, we set 
the strength of the interaction J and the Boltzmann constant fc^ equal to 1 and time is measured in Monte Carlo steps 
(MCs). After a few MCs the system enters the scaling regime in which two-time averaged quantities depend only on 
the ratio t/t w . In the inset of Fig. [3] we show the behavior of (x)c f° r different choices of t/t w . The curves are clearly 
distinct. In the main part of the figure the collapse obtained by plotting (x)c/(x) against x = (C — (C))/(f — (C)) 
is demonstrated. Moreover, the curve with t/t w = 10 is numerically indistinguishable from f{x) = 1 — x/2 for all 
x < 0. Upon decreasing t/t w one can notice a small residual dependence on t/t w , the larger the smaller t/t w . This is 
due to the fact that n is finite. By increasing n one can check that the curves converge to a master-curve behaving 
as fix) = 1 — x/2 for x < (however averaging over larger boxes reduces the extent of fluctuations and one can only 
study a small region of C values around (C)). This confirms the validity of the scaling (fT3| . 
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FIG. 2: Sketch of the interface motion in the Id Ising model. The horizontal axis represents space and the vertical one the 
Ising spin configurations. In the upper (lower) part of the figure the configuration at t w (t) are shown. 



In conclusion, for x < 0, f(x) is given by the non-interacting interface approximation. Since xq — —(C) /(l — (C)) < 
0, Eq. (fT4]) holds and the limiting can be read from the slope of f(x) in the negative x sector. The fact that 
the approximation describes very accurately the data in the full x < sector is somehow surprising since neglecting 
interface interactions is hard to justify far from x = 0. The deviation of the function f(x) from the linear shape in 
the region x > was to be expected since the approximation used to derive it, in particular the assumption 1 — C 
significantly different from zero, is not respected. In particular, since for C = 1 it must be (x)c=i = 0, /(l) = must 
hold and the curve bends downwords with respect to the non-interacting interface approximation. 

Let us stress that the non-trivial factor x((C)) in Eq. (flU)) has to be divided away to obtain X^: the slope of (x)c 
against C depends on (C), it is not constant and differs from Xoo (see the inset in Fig. [3]). Notice also that the curve 
(x)c(C, (C)) is different from x((C)) for any choice of t, t w : the two curves cross at C = (C) with a different slope. 
We shall comment on the implications of these results on time-reparametrization invariance in Sec. IIVI 



B. Ising model in d = 2, 3 quenched to T c 

We consider now the Ising model quenched to the critical temperature T c in three and two dimensions. We first 
use a Gaussian approximation and later we present numerical simulations. 



1. Gaussian approximation 
A Gaussian joint probability distribution of C and x reads 



Pg(C, Tx) = exp j-i (SC, TS X ) A (tS X )} (29) 

where SC = C — (C), Sx = X~ (x)i (C) and (x) are the mean values that we introduce as time-dependent parameters, 
and the matrix A is 

A (V C V Cx y 1 _(V x -V Cx \ 1 f3Ql 

\Vc x V x ) -{ -V Cx V C ) V C V X -VS X ■ [M) 

This PDF can only be an approximation in the critical dynamics of the finite dimensional Ising models. This forms 
is exact, instead, for quantities C and x that summed over all spins in the 0(Af) model in the large TV limit, as it has 
been analyzed in [16| . 

Deviations from the Gaussian PDF are expected for finite Af and finite N systems. Annibale and Sollich recently 
studied the critical dynamics of ferromagnetic spherical models with finite N. The detailed analysis of 1/VN correc- 
tions developed in [17| showed that these can be treated perturbatively at leading order for quenches at criticality 
and are thus amenable to analytic investigation. We shall discuss the results in this paper and how they compare to 
ours in the Discussion Section. 




FIG. 3: (Color online.) Test of the scaling hypothesis (|13p in the Glauber dynamics of the one dimensional Ising model quenched 
to T = (using N = 10 6 and I = 100, 150, 200, 125 respectively for the cases with t/t w = 2, 3, 5, 10). In the inset the quantity 
(x)c is plotted against C for different choices of t/t w given in the key and t m — 10. The heavy symbols represent the points of 
coordinates ((C), (x)) for the different choices of t/t m , which lie on the fully averaged curve x((C)) represented with a dashed 
line. The curves (x)c(C', (C)) are close to each other but there is a small visible drift for increasing t/t m and they are not 
parallel. In the main panel {x)c{C, (C)) /x((C)) , that defines the function f(x) through Eq. (|13[) . is plotted for the same set 
of data. The solid (violet) line is the linear behavior, f(x) — 1 — x/2, obtained in the non-interacting interface approximation 
of Sec. IIII A II that very accurately describes the data in the negative region of x. See the text for a discussion. 



The restricted average yields 



TfdxxP G (C,T X ) 
fd X P G (C,Tx) 



T( X ) 



V Cx C-(C) 



Vc T( X ) 



(31) 



(Note that J dC (x)c ^ (x)-) Recalling that for large t w the FDT holds T(x) = 1 — (C) in a critical quench, 
Eq. (|3ip has the scaling form (|13p proposed if the ratio Vc/Vcx is a constant. Interestingly enough, in this framework 
= —Vc/Vcx would be given by a ratio of covariances. 
The joint PDF of C and x depends on the value of £, the coarse-graining length over which these quantities are 
computed. Clearly, for I ^> Lit) fluctuations become rare and are sharply distributed around the mean values with 
Gaussian statistics. In this limit the Gaussian approximation should be very precise. However, this limit is not that 
interesting for our purposes (and for the utility of this type of measurement as a method to estimate Xoc) since the 
extent of fluctuations is heavily suppressed. In the more interesting case l~L(t) the joint PDF cannot be Gaussian 
but, still, as we shall see numerically, the approximation is of relatively good quality and yields a good estimate of 

X^. 



We have checked the quality of the Gaussian approximation and its implications on the validity of both Eq. (fie 
and Eq. (fT4|) by computing numerically Vc = (5C 2 ) and Vc x = (<5C T8x) hi the critical quench of the Ising model in 
d — 2. The angular brackets indicate here an average over the numerical distribution function. 

In Fig. [4] we display the time dependence of Vc and Vc x - The curves initially depend on time, then reach (approx- 
imately) a plateau and their ratio, shown in the upper panel, a constant. At rather long times a time dependence 
develops signaling that the Gaussian approximation becomes less accurate, a fact that was to be expected since L(t) 
increases and I is no longer larger than Lit). Focusing on the constant part of the ratio one notices that it is very 
well close to the known value of X 00l indicated with a dashed horizontal line in the figure. These findings clearly 
substantiate our conjecture in Eqs. (fT3")l and Q14p . 

Finally, let us add that the computation of V x at fixed I (where x is computed through Eq. (|2"3"|) ) yields a 
monotonously increasing function diverging in the large-t limit (not shown in Fig. [3]). We shall comment on this issue 
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FIG. 4: (Color online.) In the lower panel the covariances Vc (upper curve), and — Vc x (lower curve) are plotted against time 
for the Ising model in d = 2 quenched to T c — 2.269. The system size is C = 2 • 10 3 and the coarse-graining length is t = 85. 
In the upper panel the behavior of the ratio — Vcx/Vc is shown. The dashed line is the value of X x . 

in Sec. El 



2. Numerical results for the full model 

We start from the d — 3 case quenched to the critical point. In the inset in Fig. [5] we show the plot of (x)c versus 
C. According to the scaling in Eqs. (|15[) and (fT6|) of the fully averaged quantities in the large t w limit with fixed t/t w 
one has (C) — > and T(x) — > 1- This implies that, in this limit, the scaling (fl~3|) of the restricted averaged linear 
response is lack of content since it becomes (x)c — /(C)? just defining /. However the truly asymptotic limit cannot 
be reached numerically and, for the values of t w used in the simulations, one still finds a difference of order 0.05 — 0.15 
(depending on the curve) between (C) and 0, and T(x) and 1. Despite these important pre-asymptotic corrections in 
the independent behavior of (C) and T(x) they still satisfy the asymptotic relation (JT7J) with good precision. 

In Fig. [5] we test the scaling (|13|) for the pre-asymptotic dynamics. The collapse is excellent in the region x ~ 
and looses quality for larger values of |a;|, probably due to too important finite t w corrections. In order to check the 
conjecture l[l"3|) we evaluate the slope of f(x) in x = 0. Since x — > for (C) — ► this procedure is asymptotically 
equivalent to measuring the slope in xo , with the great advantage of using the region where the scaling (fl3|) is well 
obeyed and where we have the best statistics. The value of the slope measured in this way yields similar values for 
all the curves, all of which in remarkable agreement with Xoo ~ 0.40 — 0.43 as reported in the literature [HQ. Note 
that in the limit (C) — > the prefactor Tx((C}) tends to one and the slope of / becomes the slope of T(x)c that 
coincides with the slope of Tx((C)). The interesting regime has been compressed to a single point in the Tx((C)) 
graph but it opens up in the fluctuation analysis. The local slopes of the fluctuating relation at fixed times and the 
fully averaged relation using time as a parameter coincide. 

The situation is qualitatively similar in d = 2, although pre-asymptotic effects are stronger. This is due to the 
value of the exponent a which regulates the scaling behavior in Eqs. (fT5|) and (|16|. that, being smaller in d = 2 
than in d = 3, delays considerably the asymptotic convergence (l9j |. Indeed, in order to reduce the deviations of (C) 
and T(x) from the asymptotic values to values that are comparable to the ones found in the d — 3 case we had to 
use much longer times. Moreover, differently from the three-dimensional case, (C) and T(x) do not obey Eq. ([17)1 
pre-asymptotically as shown in the lower inset to Fig. [5] by the fact that the heavy symbols fall away from the dashed 
straight line. Due to such larger pre-asymptotic effects there is a residual time dependence in the slope of f(x) in 
the origin, as shown in the main part of Fig. [6l and the collapse of the curves is not as good in d = 3. However, 
the quality of the collapse improves as the asymptotic region is approached (namely as (C) decreases). Actually, 
the curves corresponding to the two smaller values of (C) exhibit a good collapse in the central region C ~ (C). 
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FIG. 5: (Color online.) Test of the scaling hypothesis (|13[) in the critical dynamics of the three-dimensional Ising model 
(T c — 4.511). The system size is £ = 200 and the coarse-graining length is i = 20,30,30 for the cases (t w = 20, t = 30), 
(t w = 10, t = 25) and (t w — 30, t = 50), respectively. In the inset the quantity {x)c is plotted against C. The heavy symbols 
represent the points of coordinates ((C), (x)) for the different choices of t, t w given in the key and they all lie on the FDT curve 
Tx({C)) = 1 — (C) (dashed line). In the main part of the figure the function f(x) defined through Eq. (|13[) is extracted from 
the same set of data. The solid (violet) straight line is the expected small- a; behavior f(x) = 1 — XoaX with Xoo ~ 0.41. Best 
fits of the slopes of / in x = yield 0.41, 0.42, 0.40 for the cases (t w = 20, t = 30), (t w = 10, t = 25) and (t w = 30, t = 50), 
respectively. 



Moreover, the slope slowly approaches the known value = 0.33 — 0.34 [21], which is reached in the longest run. 
We stress that also in this case the curve (x)c(C, (C)) i s different from x((C)) but, in the asymptotic limit in which 
Tx((C)) -> 1 the slope of (x)c(C, (C)) at C -> coincides with the slope of £((C)) at (C) -> and they are both 
given by -X^. 



C. Ising model in d — 2 quenched below T c 

We consider now the behavior of a ferromagnetic system quenched below the critical temperature. We restrict the 
analysis to a two-dimensional case because the task is computationally demanding and we do not expect qualitative 
differences in higher dimensions. 

We assume that a splitting analogous to Eq. (fTT))) holds also for the restricted averages, namely 

(X)c = {X)% + (X) a c 9 ■ (32) 

Since (x)c ^ s an equilibrium contribution it depends only on the time difference t — t w . Working with fixed t/t w in the 
limit t w — ► oo amounts to probe the large t limit of (x)c- Since this is a static quantity computed in an equilibrium 
state (although with a restricted average) it cannot depend on the history, and hence neither on C . Therefore, in the 
limit considered one has (x)c = X eq = /3(1 — M 2 ). 

Next, we want to show that (x)cf can be neglected with respect to (x)c; similarly to what happens for the fully 
averaged quantities. In appendix II we present a scaling argument showing that in the large t w limit (x)c = for 
a quench to T = 0. Since in this case (x)c = (x)cf an( ^ the effect of a finite temperature is not expected to change 
significantly the behavior of the aging contributions, this argument suggests that (x)cf can indeed be neglected for 
any T < T c . 

In the following we test this statement numerically. The most obvious way of computing (x)cf i s by subtracting 
(x)c f rom (x)c- However there is a by far more efficient way to compute (x)c? that consists in considering a modified 
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FIG. 6: (Color online.) Test of the scaling hypothesis (|13[1 in the critical dynamics of the two-dimensional Ising model 
(T c — 2.269). The system size is £ = 2 • 10 3 and the coarse-graining length is i = 60 except for the cases with (t w ,t) — (150, 300) 
and (t w ,t) = (50,3600) for which I = 50 and i — 85 were respectively used. In the lower inset T{x)c against C is shown 
for different choices of t and t w given in the key. The heavy symbols represent the points of coordinates ((C), (x)) and they 
fall away from the FDT curve x((C)) = 1 — (C) (dashed line). In the main part of the figure the function f(x) defined 
through Eq. (|13[) . The upper inset shows the behavior of df/dx \ x= q as C = (C) is varied. The trend seems to indicate 
limc-»o df(x)/dx — —X x ~ —0.34, as expected. 



dynamics in which flips in the bulk of domains are prevented. Since the stationary contribution is given by the reversal 
of spins well inside the domains this no-bulk-flip dynamics isolates the aging behavior with the numerical advantage 
of evolving only the small fraction of interface spins. This technique has been thoroughly tested and used in studies 
of ( X ) and (C) Hi|,[3l|. 

We have checked that also for the restricted average response the no-bulk-flip kinetics yields 
the same results as subtracting (x)c f rom (x)c- The results obtained with this kind of dynamics are shown in Fig. [7] 
One concludes that, working with a fixed t/t w (for instance, the set of data obtained with t/t w = 2 are shown in the 
figure), for any given value of C, (x)c and its slope go to zero. This guarantees that for very long times (x)c? can be 
neglected with respect to (x)c an d hence one has T(x)c(t> t w ) = 1 — M 2 and f(x) = 1, leading to a vanishing slope 
and Xoo = 0. Notice that the mechanism producing T e ff = oo in the full aging regime (and in consequence = 0) 
is the same as for the global quantities, namely the fact that the aging contribution vanishes asymptotically. 



IV. DISCUSSION 



In a series of papers it was claimed that time-reparametrization invariance is the symmetry that controls dynamic 
fluctuations [01, [TH, [H[ in aging systems with a finite effective temperature, namely low temperature glassy cases. 
A consequence of this proposal is that the fluctuating linear responses and correlation functions computed over the 
same realization (in the same subsystem of size t with the same stochastic noise) should, in the t w — > oo limit, the 
aging regime and the scaling limit 6 <C I -C C, be linked and satisfy 

X a9 (t,t w ) = x ag (C) with C = C(t,t w ). (33) 

Importantly enough, the function x~ a9 (x) is the scaling function of the global and fully averaged linear response in the 
aging regime in which the latter is a non-trivial function X as ((C)) °f the global and fully averaged correlation, see 




c 

FIG. 7: (Color online.) The quantity {x)c is plotted against C for the two-dimensional Ising model quenched to T — 1 < T c 
using t/t w =2 and different choices of t w . The system size is C = 2 • 10 3 and the coarse graining length is £ = 50, 70, 100, 150 
for the cases with t w — 50, 100, 200, 500, respectively. The heavy symbols represent the point of coordinates ((C), T(x)) for the 
different choices of t w ■ The curves are slowly turning flat as t w increases implying that Xoo — > as expected. 

Eq. JB]). The requirement of having a finite effective temperature translates into the fact that x as (C) does not vanish 
asymptotically. The statement (155]) is equivalent to saying that the fluctuations along the curve x a£, (C) are massless 
while the ones that do not follow this direction are massive and can hence be eliminated by the coarse-graining in 
the scaling limit. Time-reparamentrization invariance is a symmetry that is expected to develop asymptotically. At 
finite times £ should be scaled with a growing correlation length that diverges asymptotically. Numerical tests of 
this proposal in the low temperature dynamics of the 3d Edwards- Anderson model [131 ] , some kinetically constrained 
systems [l5| , Lennard- Jones mixtures [32[ and disordered elastic lines [33[ yielded encouraging results. 

The search for time-reparametrization symmetry in the 0{N) coarsening model in the large M limit showed that 
this symmetry is not fully realized in this case; it is, instead, reduced to global rescalings of time, t — ► Xt with A 
a constant parameter This result suggested that dynamic fluctuations in coarsening problems might follow a 

different rule although the question remained as to whether the reduction of time-reparametrization invariance to 
time-rescaling was a pathology of the large Af limit. 

Annibale and Sollich recently started the study of critical out of equilibrium dynamic fluctuations by analysing the 
ferromagnetic spherical model with finite number of spins, N, including 1 / \/N corrections 17] . The joint PDF of the 
global (summed over all spins in the system) C and T\ is Gaussian for N — > oo and the contour levels are ellipses. For 
finite N the fluctuations deviate from Gaussian statistics; however, at leading order in 1/y/N the critical fluctuations 
can be treated perturbatively, and one recovers Gaussian statistics for (C, T\) with elliptic contour levels that can be 
computed analytically. The principal axis of any of these ellipses forms an angle <fi with the C axis that is given by 

[tan^]- 1 = tan[(l/2) atan(2Vfe x /(Vb - V x ))] , (34) 

see Fig. [5] Notice that the angle <j> depends not only depends on Vc and Vc x but also on V x . As discussed by 
Annibale and Sollich, if the time-reparametrization invariance scenario applied to critical dynamics, the angle <f>, that 
is a natural measure of the slope of the cloud, should yield the fluctuation-dissipation ratio A ((C)) that relates the 
variations with time of the average susceptibility and correlation, see Eq. ([6]). In particular, it should yield when 
the long t limit is taken before the long t w limit. The analytic computation of Vc, V x and Vc x in d < 4 showed 
that the angle is not related to X in any simple way. For small time differences the variances and co-variance are 
stationary but one does not recover the FDT slope 1 from this calculation and in the opposite t 3> t w or C — > limit 
the angle tends to ir/2 meaning that the ellipses stretch in the susceptibility direction. These results invalidate one 
consequence of time-reparametrization invariance and indicate that this symmetry does not develop asymptotically 
in the critical dynamics of the ferromagnetic spherical model at leading order in 1/N. 

In this paper we analyzed the linear- response/correlation fluctuations in critical dynamics and coarsening in finite 
dimensional systems with finite dimensional order parameter. The main point of this paper is to propose the use 
of restricted averages, in which the integrated linear responses are averaged over trajectories that have the same 
value of the fluctuating two-time function C, to study fluctuations in critical and sub-critical out of equilibrium 
dynamics. We sudied the relation between restricted averaged susceptibility and fluctuation two-time function and 
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we conjectured that it yields the asymptotic effective temperature T/X^ relevant to critical dynamics. Although 
time-reparametrization invariance was not checked explicitly, our results bear some indication on the existence or not 
of such a symmetry in these cases. We summarize the results for different cases below. 

We discuss quenches to T c > first. In these case, a Gaussian approximation is rather accurate if one uses coarse- 
graining lengths, £, that are significantly larger than the growing length L(t). Clearly, as time increases L(t) goes 
beyond £ and the approximation has to be revised. The restricted average using the Gaussian PDF then naturally 
provides a slope tan 6* = —Vc x /Vc that is different from tan^>, see Eq. (|34l) . As shown in Fig. [5J this is the slope of 
the line connecting the two points of the ellipse where C takes the largest (smallest) value. This quantity does not 
depend on V x and was shown to be the one yielding Xoo . Indeed, Vc x an d Vc converge to a constant for long (but not 
too long) t and their ratio equals —X^. Since 6 yields X^ and 4> ^ 6 our results suggest that time reparametrization 
invariance does not hold in this case. Note that despite the strictly Gaussian character of fluctuations in the large- TV 
spherical ferromagnetic model, the Vc, V x and Vc x behave in a radically different way from our determination. In 
particular, the ratio between Vc x and Vc does not approach a constant in d < 4. 

V x grows very fast as a function of time when the fluctuations are computed at fixed £. Therefore, the fluctuations 
of x are very important and cannot be reduced - as compared to the ones of the corresponding correlation - by using 
a convenient choice of £. As found by Annibale and Sollich for the spherical model, the axis of the cloud tends to turn 
parallel to the x axis. This behavior is also at odds with what one would expect if time-reparametrization invariance 
were obeyed, since the cloud should asymptotically lay parallel to the function x ag (x), namely with a finite slope X^. 
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FIG. 8: (Color online.) The ellipse defining the surface of constant probability P(C, T\) = const, in the Gaussian approximation 
(|29l) . The two ellipses have Vc = 2Vc x — 1 while Vc x — 1 an d Vc% = 5 for the inner and the outer ellipse respectively. The 
slope tan 9 = — Vcx/Vc (dashed line) is the same for the two ellipses, while the slope tancj) of the major axis (continuous 
straight lines) is different in the two cases. 

The situation is different for the quench to T = in d = 1. The numerical data and the independent interface 
approximation of the kinetics of the Ising chain showed that (x)c is not equal to x(C)- Instead, we found that (x)c 
is given by a scaling function that is proportional to the non-trivial function multiplied by another non-trivial 

function of C and (C). These results suggest that time-reparametrization symmetry may hardly be realized in this 
case and that the relation between the slope of the cloud and X^ is even more hidden than in the case of quenches 
to T c > 0. Nevertheless, despite the difference between (x)c and x((C))i interestingly their relation still encodes the 
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limiting effective temperature Xoq through the scaling function / defined from (x)c/x((C))- 

For proper sub-critical coarsening the aging contribution to the restricted averaged integrated linear response 
vanishes asymptotically and, as for the fully averaged quantity T(x)c approaches the constant 1 — M 2 . In a plot 
T(x)c against C one then just sees horizontal fluctuations but this is a trivial consequence of the fact that T(x)c 
approaches a constant. This result extends the one found in [16l | for the OiAf) model in the infinite Af limit to domain 
growth with finite dimensional order parameter. 

The study of the same fluctuations could also be addressed experimentally. A recent study of the out of equilibrium 
relaxation after a quench to the Freedericksz second-order phase transition in a liquid crystal demonstrated that the 
fully averaged correlation and linear response age and are linked by an FDR with an effective temperature that is 
higher than the environmental one [34| . The analysis of the fluctuations of these quantities and the restricted average 
proposed in this paper should shed light on the generality of our scaling hypothesis and the fact that could be 
accessed by studying fluctuations. 



APPENDIX I 



Let us consider an interface / located in x(t ) (t w < t < t), namely s< = +1 for i < x(t) and Sj = —1 otherwise, 
see Fig. [5] We denote with i~(t') and i + {t') the sites surrounding / (on the left and right respectively). One has 

l-C=^, (35) 
where Ax = \x(t) — x(t w )\. Let us now evaluate the term 

1 1 N * 1 

-D(t,t w ) = -Y,Si(t) E 2 Biit>) (36) 

i=l t'—tw 

appearing in Eq. (|23|) . Let us stipulate that the velocity v(t') is +1 (—1) if the interface moves to the right (left). 
Notice that we are implicitly assuming Metropolis- like transition rates (namely the interface always moves). It is 
useful to introduce the map of accelerations at, defined as follows. Starting from <Zj = an update occurs only in the 
following cases: 

• When / starts moving at t' — t w : One changes ai-u w ) — * a i-(t m ) + 1 i a i+(t w ) ~~ * a i+(t w ) ~ 1) if the v w = v(t w ) 
is positive (negative). Here and in the following all the at not specifically mentioned remain unchanged. 

• When I stops moving at t' = t: One updates CLi+(t\ = di+(t) — 1 (aj-(t ro ) = a i-(t w ) + 1) if v (t) is positive 
(negative). 

• When the velocity changes from +1 to —1 (—1 to +1) in x(t'); One changes aj-( t n — > dj-( t M — 1, di+w) ~^ 
o»+(t') _ 1 (oi-(t') - * ai-(t') + 1) ai+(r) — * <H+{t') + !)■ 

At T = one has transition rates Wi(t') = 1 on sites i~(t'), i + {t') and Wi(t') = elsewhere. Then one has 



1 



1 for i = i-(t') 



-Bi(t') = { -1 for i = i+(t), (37) 



2 



for i£ r(t'),i+(t') 



When / moves these ±1 contributions are seeded in the region traveled which will be then summed up in D(t, t w ). If 
there are no accelerations (in the sense defined above) all these contribution sum up to zero in computing the integral 
over time in Eq. (|36[) . Since contributions to the integral come only from accelerations it is easy to prove that one has 

1 1 * 

-D(t,t w ) = E ( 38 ) 

i t'=t w 

Most of the contributions to the sum over times cancel each other. Let us indicate with tk the time when / crosses 
x(t) for the fc-th time. By definition / changes direction an odd number of times in every interval [tk, tk+i]- Since the 
contributions due to / changing direction an even number of times without crossing x(t) cancel out in the sum over 
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times of Eq. {38]), one is left only with the accelerations at (say) the last velocity reversal. Then, for each interval 
[tfc, ifc+i] there is a contribution +2 (—2) if the interface is on the right (left) of x(t) which, once multiplied by s(t) in 
Eq. (|38|) gives a contribution +2. In conclusion, for / crossing x(t) n c times one has a contribution 2n c to (l/2)D(t, t w ) 
plus the accelerations at t w and t. Since the latter yield a contribution +1 one has (l/2)D(t, t w ) = [2n c + 1 + v w ]/N 
and hence 

2 t-t w 

(D(t,t w )) c = Yl [2n c + l + v w }P c (n c ,v w ), (39) 

n c =0 «„=±1 

where Pc(n c ,v w ) is the probability of finding a particular choice of n c and v w in the restricted ensemble. 



APPENDIX II 



In this section we develop a scaling argument to compute (x)c hi a- quench to T = in d = 2. We shall assume 
Metropolis transition rates, as for d = 1 . We follow the scaling approach developed in [35| , which amounts to consider 
the relaxation of a domain of (say) down spins which at time t w has a faceted interface, as depicted in Fig. [9] The 
process ends at the final time t w + t when the domain has disappeared, namely all its spins have been reversed. The 
autocorrelation function of this process is easily evaluated as C(t w + r,t w ) = 1 — 2Nfj/N, where Njj is the initial 
number of spins in the domain. Let us evaluate the quantity (\/2)D(t w + r, t w ) in Eq. (|36[) . recalling that Si(t) = 1. 
At time t' = t w the only possible moves are the flip of a corner spin. Let us assume for simplicity that this is the 
top right, as shown in the left panel in Fig. [9] This move generates a kink which performs a random walk on the 
edge of the domain until it disappears when it reaches the boundary on the left side (or another anti-kink generated 
by the flipping of the spin on the upper left corner). In this way the first row is eliminated, and the process is then 
repeated until the domain disappears at time t' = t w + t. At T = 0, only spins the flip of which do not increase the 
energy can be updated. Hence at a generic time t' , the only contributions to (1/2)13 are those provided by the spins 
in the corners or those surrounding kinks. The contribution of the corners is always —1 while spins surrounding a 
kink contribute ±1 (—1 for the spin belonging to the domain, +1 for the other). As the dynamics proceeds, many 
of these contributions are generated, that must then be summed up in (1/2)13. In so doing, however, it is easy to 
realize that all the contributions coming from the kinks sum up to zero, since they always occur in pairs. Hence one 
is left with the contributions from the corners only. As shown in the right panel in Fig. [5] topological reasons fix the 
number of corners to be 4 at all times. The contribution —1 of these spins lasts for all the time r of the process. Then 
one has (1/2)13 = —At/N. The next step is to evaluate (D)c- Since C is univocally determined by No, computing 
(D)c simply amounts to determine the average time needed for a faceted domain of Nd spins to disappear with 
zero temperature dynamics, namely (D)c = —8(t)n d /N. In [35J] it was shown that (t)n d = No/A. Then one has 
(D) c = -2N D /N = -(1-C), and from Eq. {23]): 

(X)c(t,t w )=0. (40) 
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FIG. 9: The relaxation of a domain with faceted interfaces at T — 0. Left: The initial configuration, where only corner spins 
can be flipped. Right: Configuration at a generic time t 1 with a kink. The four corner spins producing the relevant contribution 
to D are evidenced. 
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